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Abstract 

The decay of (disorder-averaged) static spin correlation functions at T = 
for the one-dimensional spin-1/2 XX Z antiferromagnet with uniform longi- 
tudinal coupling J A and random transverse coupling JA, is investigated by 
numerical calculations for ensembles of finite chains. At A = (XX model) 
the calculation is based on the Jordan- Wigner mapping to free lattice fermions 
for chains with up to N = 100 sites. At A / Lanczos diagonalizations are 
carried out for chains with up to N = 22 sites. The longitudinal correlation 
function (SqS*) is found to exhibit a power-law decay with an exponent that 
varies with A and, for nonzero A, also with the width of the Aj-distribution. 
The results for the transverse correlation function (SqS^:) show a crossover 
from power-law decay to exponential decay as the exchange disorder is turned 
on. 
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The combination of randomness and quantum fluctuations is well-known to be a fer- 
tile ground for interesting physical phenomena including Anderson localization. In one- 
dimensional (ID) tight-binding systems the rule is that disorder always leads to localiza- 
tion. However, if the randomness is purely off-diagonal, the localization length diverges at 
the band center,0 which is bound to affect the decay law of correlation functions. One par- 
ticular off-diagonally disordered fermion model, the ID half-filled tight-binding model with 
random hopping, is equivalent to the special XX case (A = 0) of the ID s = 1/2 XX Z 
model with random transverse exchange coupling, described by the Hamiltonian 

H = JXJMW! + S?S? +1 ) + AS*S* +1 ] . (1) 

i 

The uniform longitudinal spin coupling corresponds to a fermion interaction. Here we con- 
sider the range < A < 1, use periodic boundary conditions, and take the random transverse 



coupling JXi to be described by a Gaussian distribution with Aj = 1, AjAj = 1 + cr 2 (%. The 
spin correlations at T = of this model were recently investigated by means of a real-space 
renormalization group (RSRG) methodic based in part on ideas from earlier work,@ and by 
means of a finite-chain study! 

One interesting proposition made in the context of the RSRG study is the existence of 
a random- singlet phase with algebraically decaying spin pair correlations :!'! 

(S«S?) ~ (-l) p r-*», a = x,z. (2) 

The singlet nature of that phase would imply that the characteristic exponent r] a assumes 
the same value in the longitudinal (z) and transverse (x) correlation functions, in marked 
contrast to the case with no exchange disorder (<r — 0), for which we know the exact result!! 

rj x = l/r] z = 1 — arccos A . (3) 

The RSRG study further predicts that this exponent value is r\ x — rj z — 2, independent of the 
longitudinal coupling A and the disorder strength a, provided the latter is not too small. 
It is indeed quite unusual that the anisotropic randomization of an anisotropic exchange 
interaction should effectively remove the effects of anisotropy in the spin correlations. 



Here we report results of a finite-chain study which goes significantly beyond that of Ref. 
[5] in statistics and system sizes. For the XX model (A = 0), we carry out the computation 
in the (free-) fermion representation, which enables us to handle chains with up to iV = 100 
spins and beyond. For A ^ we must resort to Lanczos diagonalizations. Here the largest 
system for which we can perform the computation with reasonable statistics has N = 22 
sites. For graphical purposes we shall consider, henceforth, the absolute value, \(SqS^)\, of 
the spin pair correlations. 

We first consider the case A = (XX model). If the longitudinal correlation function 
does exhibit power-law decay, \(SqS*)\ ~ r~ Vz , as predicted, then the exponent r\ z also 
governs the iV-dependence of the function \{SqS z n , 2 )\ in a cyclic chain of N sites.H We have 
evaluated this quantity for systems with N < 100 sites and for disorder strengths a < 2, all 
with ensemble averages over up to 10 5 configurations. 

The data analysis yields r] z = 2 independent of a . This is consistent with the RSRG 
prediction^! but in contradiction to the earlier finite-size studyi where a significant o- 
dependence of r] z was observed.B Our data also confirm that the disorder-averaged logarithm 
of \(SqS*) \ exhibits the decay law ~ exp(— r 1 / 2 ) as predicted in Ref. ^. 

The decay of the transverse correlation function (SqS^) is much more sensitive to the 
presence of exchange disorder, as is demonstrated by the data shown in Figs. [I] and |2|. 
In the main diagram of Fig. |l] we show the function \(SqS*)\ versus r in a logarithmic 
plot for ensembles with different disorder strengths. Turning on the exchange disorder with 
gradually increasing a causes the transverse correlations to decay more and more rapidly as 
one might expect .! 

For < o < 0.5 the data describe a power-law behavior. This is also evident in the 
bundle of curves near the top of Fig. [| which shows the r-dependence of the function 
| (5*0 Sf) | semi-logarithmically at o = 0.4 for various system sizes. The values \{SqS^^ 2 )\ at 
the endpoints of these curves plotted vs N/2 in a log- log graph fall onto a straight line, and 
the slope of that line determines the exponent rj x . This is illustrated by the full squares in 
the inset to Fig. |[ 



The a-dependence of r\ x as obtained from this procedure is shown in the inset to Fig. [T|. 
For the system without exchange disorder we reproduce the exactly known value r] x = 1/2,1 
which is a special case of ([|). As a increases from zero, r] x grows gradually and monotonically, 
at first slowly, then more and more rapidly. 

For a > 0.5 the curves in Fig. [I] suggest the occurrence of a crossover from algebraic 
decay to exponential decay, \{SqS*)\ ~ exp(— r/£), in the range of r for which we have data. 
The exponential character of the decay is more strikingly manifest in the lower bundle of 
data shown in the main plot of Fig. |2|, representing the function | (SqS^:) | at a — 1 for various 
system sizes.! The smallest expectation values are known only with considerable (relative) 
uncertainty despite the augmented statistics. 

The triangles, which represent the values |(SoS^y 2 )l vs N/2 in this semi-logarithmic plot, 
are consistent with a straight line. Its slope determines the disorder-induced correlation 
length £. Over the range of disorder strengths, where our data suggest exponential decay of 
| (So 5^) | , £ thus determined decreases monotonically with increasing a. 

Our data for the exchange disordered XX model are consistent with two alternative 
scenarios, which are equally interesting: (i) There exists a transition at some nonzero value 
of the disorder strength, a c ~ 0.5, from algebraically to exponentially decaying transverse 
spin correlations, (ii) A transition of the same nature occurs at o c = instead, which 
produces very similar crossover effects in the finite-chain data. A more extensive study for 
longer chains and with better statistics will be necessary to discriminate with confidence 
between the two scenarios.0 The data are definitely incompatible with a persistent power- 
law decay as predicted by RSRG. 

Now we turn to one case, A = 0.75, with fermion interaction (XXZ model). Since the 
computations are much more involved, the available data are limited by comparison with 
the case A = 0. At A ^ neither our data for the longitudinal correlations nor those for 
the transverse correlations are compatible with the RSRG predictions. 

The function |(SqS*)| for various system sizes and a = 1.5 is shown logarithmically in 
Fig. |||. The endpoint data (r = N/2), which fall neatly onto a straight line, describe a 
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power-law decay with rj z = 1.26. The exponent values obtained for two smaller disorder 
strengths are r\ z = 1.17 (a = 0.5) and r\ z = 1.31 (a = 0.25). The exact result (|3|) for a = 
assumes the value r\ z = 1.298 . . . 

All combined, the data suggest that the function \(SqS*)\ is governed by a power-law 
which persists in the presence of randomness. The cr-dependence of the exponent rj z appears 
to go through a minimum of considerable depth at o ^ 0,0 which implies the curious 
phenomenon that the longitudinal correlations are enhanced by a small amount of transverse 
exchange disorder relative to the correlations in the uniform-exchange system. 

The data for the transverse correlations \(SqS^)\ at A = 0.75 exhibit properties very 
similar to what we have observed and described for the free-fermion case (A = 0). For not 
too large disorder strengths (a < 0.5), we see a power-law behavior with an exponent that 
increases monotonically from the exactly known value i] x = 0.769 ... at cr = 0, as given by 
expression (|3|), to r] x = 1.00 at a = 0.25 and rj x = 1.49 at a = 0.5, at which point a crossover 
to exponential behavior makes itself felt. The exponential decay law at a = 1.5 is quite 
evident in the semi-logarithmic plot of Fig. [|. 

The discrepancies between our results and the RSRG predictions of Refs. |2],|3| call for an 
explanation in future studies. Possibly, the strongly anisotropic nature of the exchange in 
the model system ([]]) - even for A = 1 - is not adequately taken into account by the RSRG 
procedure, which derives from a method originally developed for a model with isotropic 
exchange E 

In order to gain further insight into the properties of the XXZ chain with random 
exchange, we plan to investigate the nature of low-lying excitations and the properties of 
dynamic correlation functions. At A = the Jordan- Wigner mapping to free fermions will 
make it possible to carry out these calculations for large systems. At A / the KPM 
method developed recentlyB promises to be an adequate calculation! instrument. 

The work at URI was supported by NSF Grant DMR-93-12252 and by the NCSA at 
Urbana-Champaign. 
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FIGURES 

FIG. 1. Log-log plot of \(S$S*)\ at T = in the random-exchange XX chain (A = 0) with 
N = 40 spins at disorder strengths a = 0, 0.1, . . . , 1 (top to bottom). For < a < 0.5 the 
expectation value of (SqS*) has been averaged over 10 4 configurations, and for a > 0.5, over 10 5 
configurations. The inset shows the (effective) decay exponent function of the disorder 

strength a. 

FIG. 2. Semi-logarithmic plot of \(SqS^}\ with r < N/2 at T = in the random-exchange XX 
chain (A = 0) with N = 18, 22, 26, 30, 34, 40, for a = 0.4 (upper set of curves, averaged over 10 4 
configurations) and a = 1 (lower set of curves, averaged over 10 5 configurations). The data points 
for a = 1 at maximum distance (r = N/2) are marked by full triangles. The straight line which 
best fits these data points is shown dot-dashed and determines the correlation length £. The inset 
shows the data for a = 0.4 in a log-log plot. The data points at maximum distance (r = N/2) are 
marked by full squares. The straight line which best fits these data points is shown dot-dashed 
and determines the correlation exponent r\ x . 

FIG. 3. Log-log plot of \(S§S*)\ for the XX Z model with uniform longitudinal exchange 
(A = 0.75) and random transverse exchange (a = 1.5) on chains of various lengths. The data 
points represent averages over 1000 configurations for N = 6, ... , 18, and 450 configurations for 
N = 22. 

FIG. 4. Semi-logarithmic plot of \(SqS^)\ for the XX Z model with uniform longitudinal ex- 
change (A = 0.75) and random transverse exchange (a = 1.5) on chains of various lengths. The 
data points represent averages over 1000 configurations for N = 6, . . . , 18, and 450 configurations 
for N = 22. 
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